arXiv: 1503.05387vl [math.CA] 3 Feb 2015 


Asymptotic properties of biorthogonal polynomials 
systems related to Hermite and Laguerre polynomials 


Yan Xu 

School of Mathematics and Quantitative Economics, Center for Econometric analysis and 
Forecasting, Dongbei University of Finance and Economics, Liaoning, 116025, PR China 


Abstract 

In this paper, the structures to a family of biorthogonal polynomials that ap¬ 
proximate to the Hermite and Generalized Laguerre polynomials are discussed 
respectively. Therefore, the asymptotic relation between several orthogonal 
polynomials and combinatorial polynomials are derived from the systems, which 
in turn verify the Askey scheme of hypergeometric orthogonal polynomials. As 
the applications of these properties, the asymptotic representations of the gen¬ 
eralized Buchholz, Laguerre, Ultraspherical (Gegenbauer), Bernoulli, Euler, 
Meixner and Meixner-Pllaczekare polynomials are derived from the theorems di¬ 
rectly. The relationship between Bernoulli and Euler polynomials are shown as 
a special case of the characterization theorem of the Appell sequence generated 
by a scaling functions. 
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1. Introduction 


The Hermite polynomials follow from the generating function 


= E 




z G C, x G ' 


m =0 


m\ 


which gives the Cauchy-type integral 
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H m (x) = — & e xz ~~z-^+^dz. 
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( 1 . 2 ) 


The derivatives of the Gaussian function, G(x) = 






, produce the Her¬ 


mite polynomials by the relation, (—l) m G( m )(:r) = H m (x)G(x),m = 0,1,_ 

Therefore the orthonormal property of the Hermite polynomials, 


1 

m\ 


H m (x)H n (x)G(x)dx = 6„ 


can be considered as a biorthogonal relation between the derivatives of the 
Gaussian function, {(—l)"G < - n l : n = 0,1,...} and the Hermite polynomials, 
: m = 0 , 1 ,...}. 

The Hermite polynomials have been extensively studied since the pioneer 
article of C. Hermite Q in 1864. It has many interesting properties and appli¬ 
cations in several branches of mathematics, physical sciences and engineering. 
A rich source of orthogonal polynomials, for instance, Gegenbauer Q, La- 
guerre [3|, Charher flfll. Jacobi [l^, Meixner-Pollaczek, Meixner, Krawtchouk 
and Hahn-type polynomials Q have asymptotic approximations in terms of 
Hermite polynomials, which is known as the famouse Askey schemej2ll l25l |20|. 
The asymptotic relations among other hyp ergeometric orthogonal polynomials 


and its q-analogue can be found in 


25[. The asymptotic representations of 


other families of polynomials, such as the generalized Bernoulli, Euler, Bessel 
and Buchholz polynomials are also considered B flbl. I 26 1. 

In 2?||, S. L. Lee extended the biorthogonal properties between the deriva¬ 
tives of the Gaussian function and Hermite polynomials to a family of scaling 












functions with compact support and a family of Appell sequences which ap¬ 
proximate to the Gaussian function and Hermite polynomials respectively. The 
Appell polynomials are also called scaling biorthogonal polynomials which are 
eigenfunctions of a linear operator and the distributional derivatives o f <j> are 
the eigenfunctions of its adjoint corresponding to the same eigenvalues 12N . In 
particular, the Appell polynomials generated by the uniform B-spline are the 
classical Bernoulli polynomials which asymptotic approximate to the Hermite 
polynomials by suitably normalized 271. 


The main objectives of this paper is to extend these properties to a family of 
non-scaling functions that approximate to the generating functions and to con¬ 
struct a family of biorthogonal polynomials that approximate to the Hermite 
polynomials and Laguerre polynomials respectively. The asymptotic properties 
between Hermite, Laguerre and other orthogonal polynomials which are known 
as Askey Scheme are derived from these theorems as simple cases. The rela¬ 
tionship between the Appell sequence polynomials and the a-scaling compact 
support functions are also considered. 

This paper is organized as follows: In section 2, we present the framework 
of the biorthogonal polynomials system which related to Hermite polynomials 
and Gaussian function. The characterization theorem of the Appell sequence 
generated by the scaling functions are also shown in this section. In section 3, as 
the applications of senction 2, the asymptotic relations among Bernoulli polyno¬ 
mials, Euler polynomials and B-splines are studied. The new identical relation 
between Bernoulli and Euler polynomials are shown as a special case. In sec¬ 
tion 4, we generalize the asymptotic relationship to a family of hypergeometric 
orthogonal polynomials related to Hermite polynomials. The asymptotic prop¬ 
erties of Generalized Buchholz polynomials and Ultraspherical (Gegenbauer) 
polynomials are considered as the applications. In section 5, we generalize the 
biorthogonal systems to the Laguerre polynomials and derive several asymptotic 
properties in Askey scheme. 
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2. Biorthogonal polynomials approximate to Hermite polynomial 


Let C°°(R) denoted for the space of infinitely differentiable functions. If 
4> : C'°°(R) —> R is a linear functional, we shall write (<f>, v) = v G C°°(M). 
The linear functional <j> is continuous if and only if there is a compact subset K 
of R, a constant C > 0 and an integer k > 0 such that 


|| < Cmax sup \v^(x)\. 
i<k xGK 

We denote the space of distributions with compact support by £'(R). Integrable 
functions and measures with compact supports belong to £'(R). If / is a com- 
55 pactly supported integrable function then it is associated with the distribution, 
which we still denote by /, defined by 

(/» : = [ v(x)f(x)dx, v G C°°(R). 

J R 

If m is a compactly supported measure on R, then it is associated with the 
distribution, which we still denote by m, defined by 

(m, v) \= j i/(x)dm(x), v G C'°°(R). 

Jr 

Any (j> G £'(R) has derivatives </>(") of any order n and they are defined by 





n = 0 , 1 ,.... 


Taking a compactly supported distribution <fi G £'(R), then for any integer 


n > 0 , 

= (-1 r (<j>,z n eU z ) = (-1 ) n z n t(iz), 

6 o where <j>(-) denote the Fourier transform of <j>{x). 

if m ^ o, 

(Vi r^\£^-\ = z n (2.i) 

in a neighborhood of 0. Since <j> is compactly supported, (f> is analytic. So we 
can define a sequence of polynomials, Pm, by the generating function 
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XZ 




E 

m —0 


Pm{x) 

m\ 


( 2 . 2 ) 


4 



(2.3) 


It follows from EH and (EH that for any integer n > 0, 

zn = jr A-i) n 0 (n) , ^M\ 

m=0 ' m. / 


which gives the biorthogonal relation 


=5 

ml / 


(2-4) 


Definition 2.1. A sequence of polynomials, {P m {x) : m € N}, is an Appell 
sequence if P m (x) is a polynomial of degree m and 

P'mW = mP m- 1 - 

Differentiating (12.21) with respect to x and equating coefficients of z m in the 
resulting equation gives 


P! rn (x)=mP rn - i(x), m=l,2,..., 


(2.5) 


which implies P m (x ) are Appell sequence of polynomials. Therefore the distri¬ 
bution (j) generates an Appell sequence of polynomials by the generating function 


We consider a family of sequences of biorthogonal polynomials, { Pn,™ '■ rn = 
0,1,..., N = 1.2 ...}, that are generated by a sequence of functions, <pN , which 
converges to the Gaussian function 




z _ Pn, m(x) 

^ ml 


( 2 . 6 ) 


m—0 


Definition 2.2. Let (f>N be the standardized form of f>N, 

4>n(x) = a N f) N (a N x + pn), 

where pn and a jy are the mean and variance of (f>N- Define the standardized 
form of the biorthonormal polynomials, Pn,™, of {Pn,™ : to = 0 ,1 ,...} by 


Pn,™{x) = <J N m PN,m(&NX + pn)- 


(2.7) 


Then the following biorthogonal relations for the standardized biorthogonal poly¬ 
nomials follow from 

' (- 1 ) n ^7v ) > Pn,™) = S m , n ym,n> 0. 
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Further, the generating functions of PN.m are given by[27 


PN,m{x) 


( 2 . 8 ) 


4> N (iz) ~o 

Theorem 2.1. Let 4>n{x) satisfy the following conditions: 

70 (1) There exist constant r > 0 such that for any e, there is a sufficient large 

Nq, for any N > No, it holds 


4> N ( iz ) - e " 2 


z < r. 


(2.9) 


(2) Let {PN,m(x) : m = 0,1,...} he the biorthogonal polynomials generated 
by the functions, <f>N, as in 12.81) . 

Then for each m = 0,1,..., PN,m(x) converges locally uniformly to the Her- 
75 mit polynomial H m (x) as N goes to infinity. 

Proof. Since cj) N ( 0) = 1, we can choose a neighborhood U of the origin so that 


4>(iz) > ^ and e 2 
at 0 and radius r so that (ED is satisfied. 
Noting 


> tj for all z € U. Take a circle C inside U with center 


e xz _ v- P N , m {x) ~ H m {x) „ 
£ — Z 


(j> N {iz) e 2 m=0 


TO! 


( 2 . 10 ) 


The coefficients of the Taylor series (12.101) are represented by the Cauchy’s 
integral formula 

6 f „\ zj r„\ _ m! I e xz (e*r - 4> N (iz)) 

Pn, m\%) Hm\pC) — _ . V ^ „ dz. 


2n i 


c z m+1 (j) N {iz)e z 2 


Noting 4>n(x) satisfy the condition (1), which means 3 r > 0, for a real number 
A > 0, there is a sufficient large No, for any N > No, it holds 

A 


<t>N ( iz ) - e 2 


< 


— ? 


\z\ < r. 


CTN 


( 2 . 11 ) 
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Therefore 


PN,m{x) Pm{x] 


< 


< 


< 


to! 

r |e xz | 

e‘ 2 -<t> N (iz) 

27r J 

to! 

C j.m+1 

( 6 

4>ivfa) 

xRe(z) A 

O N 

z 2 

e~ 

2tt J 

4(to 

C j-m+1 

)e rx A 

4>iv( iz ) 

*2 
e 2 


I dz\ 

\dz\ 


cr nt 


Since <tat —> oo as N —> oo, It follows that for each to, PN,m(x) H m {x) 
uniformly on compact sets. □ 

2.1. Appell sequence generated by scaling functions 
as The refinement equation 

(j>n{x) = / a4> n {ax - y)dm n (y), x G R, n = 1,2,..., (2.12) 

J R 

where a > 1 and {m n } is a sequence of probability measures with finite first 
and second moments. Equivalently, (12.121) can be expressed in term of Fourier 
transforms in the frequency domain in the form 


$n{n) = fhn , h 


(2.13) 


Theorem 2.2. If two Appell sequence polynomials, P m (x) and Q m (x ) are gen¬ 
eralized by 

p xz 00 ~ra 

y— =E P ™(^) —, ( 2 - 14 ) 

- rn\ 


m —0 


and 


p xz _ _ m 

y— = E Qm{x )~T 
<j>(iz) - m\ 


(2.15) 


m =0 


respectively. Then <j)(x) is a a-scaling compact supported function with mask 
tp{x) if and only if 


E 

k =0 


(X 1 ^ J Pk (ox) Qrn—k (ox) Q m (2x). 


(2.16) 
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Proof. Since 


41 J2^- m Pm(x) Z -r, ^-=J2a-™Q m (x) Z -, 


' o: ' m =0 


then 




— [E a m P m { x ) m\ ) ( E a m Q m ( x ) m\ I ( 2 - 17 ) 


\j =0 

00 / m 


\m —0 


ml 


= E E“" m 


1 ! ■ 


( 2 . 18 ) 


ra—0 \/c—0 

Suppose that the Appell sequence polynomials, P m {x) and Q m (x), satisfy (12.161) . 


E< 

k—0 


m 


Pk {olx) Qrn—k (^3?) — Qm(2x). 


Then 

00 / m 

EE- 

m—0 \fc—0 
90 Therefore 


At-, f ax \ ^ /ax\\ z m , 

7 / m=0 


^(f)<Xf) 


00 / m 

= EE 

m=0 \k—0 


a 


= E = 

Z - 7 7T7.I 


m—0 


w! <j,{iz) 


m\ ^ /ax\ _ /ax\ \ z . 

Pfc (t) Qm ~ k (t)J ( 0) 

( 2 . 21 ) 


We see that 


${—)H~) = 

a a 


which imply <j>(x) is scaling function with mask ip(x). 

Suppose that </>(x) is a a-scaling compact supported function with mask 
ip(x), satisfying the scaling equation (12.121) . Then the Fourier transform <fi is 
given in (12.131) shows that 

^ ( a") ^ ( ck) = ( 2 - 22 ) 
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It follows that 

OO 

e Qm(2x) 


„2xz 


„2xz 


m =0 


m! 


Therefore 


j(iz) £(£)£(£ 

/ OO 

= E«- 


Vm=0 


~ra \ / 00 ~ra \ 

Pm (ax) — E a~ m Qm ( ax) — 

m 7 v^o TO 7 

™ /to\ „ , N „ . A 2™ 


00 / m / \ \ m 

= E I E a ( k ) ^ aX ^ ^ m_fc ( a;c ) ) ^TT' 

m—0 \k—0 ^ ' / 

m / \ 

E a_m ( >, ) p k ( a 7 Qm-k (ax) = Q m (2x). 

fc=0 ' ' ' 


□ 


3. Generalized Bernoulli polynomials, Euler Polynomials and B-splines 

B-splines with order N, which is denoted as Bn(-), is defined by the induc¬ 
tion as 

( 

1 if a; € [0,1), 


Bx(x) = 


0 otherwise 


and for IV > 1 


Bpj = B\ * B 


n- 1, 


where * denotes the operation of convolution which is defined by 

(• + 00 


/ -t -00 

f(t - y)g(y)dy , 

-00 


for / and g in L 2 (R). 

The Fourier transform of B^(x) is 


Bn(<x) = 


1 — e 


1 U) 


N 




Bj\f (x) also satisfies the scaling function as follow: 


N 


Bn(x) = 2 E ( N ) b n(2x - j), 

3=0 Z 


(3.1) 


9 



where the mask ipit(k) := ^( k )- Equivalently, (13.11) can be expressed in term 
of Fourier transforms in the frequency domain in the form: 


B n {w) = ^jv(w/2)Bat(u;/2), 

- / -i^\ N 

wo where the Fourier transform of the mask is ipNiui) = I — J 

The asymptotic properties of B-splines have a long history going back to the 
physicist Arnold Sommerfeld who showed that Gaussian function can be approx- 

In 1992, Unser and his colleagues 13 


Q. 


imated by B-splines point-wise in 1904 | 
proved that the sequence of normalized and scaled B-splines tends to Gaussian 
105 function in L p space as the order N increases. L. H. Y. Chen, T. N. T. Goodman 
and S. L. Lee[l 6 | considered the convergence orders of scaling functions which 
asymptotic to normality. A result due to Ralph Brinks H generalized Unser’s 
result to the derivatives of the B-splines. Yan Xu and R. H. Wang [29J gave the 
convergence orders of the approximation processes and showed the asymptotic 
no relationship among B-splines, Eulerian numbers and Hermite polynomials. 


Theorem 3.1. '2t] Let be k £ N, for N > k + 2, the sequence of the k-th 
(k) 

derivatives, B N , of the B-spline converges to the k-th derivative of the Gaussian 
function 

fc + 1 / /- \ 

1 


D K exp - 


and 



O 


(3.2) 


(- 1 )* 

V2n 


H k (x)G(x), (3.3) 


where the limit may be taken point-wise or in L p (M),p £ [ 2 ,oo). 

iflE 


Generalized Bernoulli 


"1 


10j and Euler polynomials [22, m of degree m, 
order N and complex argument 2 , denoted respectively by B^(z) and E^(z) 


can be defined by their generating functions, 


(e“ - 1 )^ 

2 N e uz 
(e“ + 1 )^ 


_^B»(z) .. 

- 2^ — 


m =0 
oo 


m! 

7 n, 


_ ^2 


m —0 


ml 


M < 27T, 


M < 7T. 


(3.4) 

(3.5) 
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115 In paper 


27[, S. L. Lee has proved that the Appell polynomials generated 


by the uniform B-spline of order N are the generalized Bernoulli polynomials of 
order N and when suitably normalized they converge to the Hermit polynomials 
as N —> oo. Since the B-splines approximate the Gaussian functionlld Il7l 124 


27 


29(, they can also be used as a filter in place of the Gaussian filter for linear 


120 scale-space 


3C . 


Corollary 3.1. '27] 



N—>oo 




Proof. Recall that the uniform B-spline, Bn{x), of order N, is the scaling func¬ 
tion satisfying 

N 

B n (x) =2^T-^ w (jB N (2x-j), 

j=0 ' J ' 

where the mask i]>N(k) := ^-()(). Equivalently, the scaling equation can be 
expressed in term of Fourier transforms in the frequency domain in the form: 


Bn(u) = i/; n {w/2)B n {w/2), 

- /. 

where the Fourier transform of the mask is —) • From the 

fourier transform of Bn we have 

( e “ - l\ N 

B N (ioj) = f——— J . (3.6) 

By (12.21) . (12.41) and the generating function of Bff(z), the generalized Bernoulli 
polynomials, {Bff(z) : m = 0, 1,...}, are biorthogonal to the derivatives of the 
B-splines , B^\ 

((-1 ) n B^(z), = 6 m , n . (3.7) 

The standardized B-splines, 


B n (x) = 


B 


N 


N N 
12 X+ 12 


converges uniformly to the Gaussian function, G{x), and an estimate of the rate 
of convergence is given in 


m 


27 


29| . 
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By Theorem 12.11 we have 


125 


lim 

N—>oo 




H rn i K z'). 


□ 


It is well known that the binomial distributions converge to the normal 
distribution in the sense that 



where Xn = VNx/2 + N/2. Let < 7 jv = \fNx/2 and /jn = AT/ 2 , then the 
standardized binomial distributions, := o’n4’n(<xnX + hn), converges 

uniformly to the Gaussian function G(x) = -^== e~ x / 2 . 

The Appell polynomials generated by the binomial distributions are the 
130 generalized Euler polynomials of order N and when suitably normalized they 
converge to the Hermit polynomials as N — > oo. 

Corollary 3.2. 

/ 4 \ 9 Jn N 

^ U) ) = W (3 - 9) 

Proof. From the Fourier transform of ^n{x) we have 

/s / -U t \ N 

■ (3-10) 

The Appell polynomials generated by binomial distributions, ^ n (fc) := 
are the generalized Euler polynomials, E^{z), m = 0,1,..., that are biorthog- 
onal to the derivatives of , iI>n{z), 

((-1 ) n 1>P(z), = 5 m ,n- (3.11) 

By Theorem (ED, the normalized generalized Euler polynomials E^(z) con¬ 
verge to the Hermit polynomials, H m [x), as IV —> oo, 

/4\? „ r JJj N 

JTooU) + = (3 - 12) 

□ 


12 





Corollary 3.3. Generalized Euler and Bernoulli polynomials satisfy 


BZ{z) = ^Y J ( V f] E ^{z)B^_ k {z). (3.13) 

k =0 ' ' 

Proof. Let 4>n(w) = — J and B n (lu) = — J • By (13.41) (13.51) . we 

can see that the Euler and Bernoulli polynomials are generated by and 

135 B N {iu) respectively. Recall that the uniform B-spline, Bjv, of order N, is the 
scaling function satisfying 

Bn{<*>) = iPn{oj/2)B n (uj/2), 


and the Fourier transform of Bn is 


B n (uj) = 


1 — e 


ICO 


N 


U! £ R. 


Therefore by Theorem 12.21 it holds 

k—0 ' ' 


(3.14) 

□ 


4. Generalized Buchholz polynomials and Ultraspherical (Gegenbauer) 
ho polynomials 


Theorem 4.1. Let {PN,m{x) ■ rn = 0,1,...} be the polynomials sequence 


generated by o 


OO Pn ,m{x)z' n 


There are constants r > 0 and 


A > 0, such that for all sufficient large N, it holds 


and 


In( 


x, z — e 


< —, 
— rr t 1 


(p N {i z ) — e *2 < -A-, for |^| < r. Then for each m = 0,1,..., PN,m{x) con¬ 

us verges locally uniformly to the Hermit polynomial H m {x ) as N goes to infinity. 


Remark 4.1. When fjq{x,z) = e xz , theorem (£ZP tuims to theorem m- 


Proof. Since 4>n( 0) = 
f(iz) 


> \ and 


e 2 


1, we can choose a neighborhood U of the origin so that 
> i for all z £ U. Take a circle C inside U with center at 


13 















0 and radius r, so that for all sufficient large N, it holds /n{x,z) — e xz < 

'Z z 2 A 

150 and (j> N {iz ) — e ~ < ^-, for \z\ < r. 

The coefficients of the Taylor series 

f N (x,z ) e xz ^ P N ,m(x) ~ H m (x) m 

— -p- = y. -j- 2 i 4 - 1 ) 

<Kiz) e T m! 

are represented by the Cauchy’s integral formula 

- , . tt t \ rn! / iffisr(x,z) e xz \ 

P N,m{x) - H m (x) =-i ^5—--7 ^ 

j7r * J c z y ^(iz) e = / 

TO ! / (e^/jv(a:, 2) - e xz j>(izf) 


2m Jr z m+1 


c z m+1 (j)(iz)e 2 


Therefore, for any given real numbers A > 0, there is a sufficient large jVq, 
for any N > Nq, it holds 


PN,m{x) H r n(x') ^ 


! r e a 2 f N (x,z)-e xz 4>(iz) 


c \ 7 m +II 


1 r e~i f N (x, z) — e xz e~i + e xz e ~2 — e xz f>(iz) 


r m+ 1 e ~2 


m l r e Re ^ f N {x,z)-e xz + e xRez - <j>(iz) 

< w 1 f -~- \ dz \ 

Jc r m+1 (j)({z) I 


A{m\)A 


Since <jn —> 00 as N — >00, It follows that for each to, pN,m(x) _■* H m {x) 
155 uniformly on compact sets. □ 


Generalized Buchholz and Ultrasplrerical (Gegenbauer) polynomials of de¬ 
gree to, order N and complex argument x, denoted respectively by Pjj ( x ) and 
Cj[(x), can be defined by their generating functions, 


O x(cotz —^)/2 


1 wo / sin Z 




14 



and 


(1 - 2xz + z 2 )~ N = Y, CZWz m , -1<x<1, M < 1. 

n =0 


leo Buchholz polynomials are used for the representation of the Whittaker functions 
as convergent series expansions of Besse functions [fj . They appear also in the 
convergent expansions of the Whittaker functions in ascending powers of their 
order and in the asymptotic expansions of the Whittaker functions in descending 


powers of their order 


i65 be found in [l 1]. 

There are well known limits 


12| . Explicit formulas for obtaining these polynomials may 


13 


lim (A) 2 P%(-2VWx) = —H m (x) 
N—too \ N J m\ 


and 


lim (2 N)~^C» 
iV—>oo 


d = 


1 




.V2NJ rnV 

These limits give insight in the location of the zeros for large values of the limit 
parameter, and the asymptotic relation with the Hermite polynomials if the 
parameter N become large and x is properly scaled. 

Many methods are available to prove these and other limits 15]. We can get 
these asymptotic results from theorem 14. II as simple cases. 

Lemma 4.1. For any \z\ < n, 
lim sinc N 

Proof. Set 


N—>oo 



Ln{z) '■= fVTn 


Then with zn = it holds 


smc 



(4.2) 


Ln{z) := IV In 



In 


= 3z 2 


sinc( 


£ / 12 ' 
2 V JV' 


3 z 2 /N 


— = 3z 


2 In [sinc( 2 jv)] 


V N 


15 













Since it holds sinc(O) = l,sinc^(0) = 0 and sinc^(O) = —for the limjv^oo, 
175 and hence zn —> 0, we may apply L’ Hopital’s rule twice: For any \z\ < n, we 
have 


lim Ln(z) 
N—>oc 


sine ^(zn) 


O 2 r sinc (1) (zjv) 0 2 

6z lim -;—-r- = 6z hm - , , 

n-hx> 2zNSmc(ZN) N—too 2 sinc( 2 jv) + 2zjvsinc ( ' i - ) (zn) 


3z- 


sinc^ (0) 


z 

Y 


It follows: For any \z\ < it, 


lim sinc^ 
N—> oo 




□ 


Corollary 4.1. 


lim , 

N—>oo V N 


P£(-2Y3Nx) = —H m (x). 

ml 


(4.3) 


Proof. Let ajy = \/ then 


oo / \ m / 

E (4) =e- 6 ' /5l ™("^ 

m— 0 \V2/ 


V 2 ct N \ ( sin 


N 


V2(tn 


Let f N (w,x) = e - 6v5ia,, ( cot *-^) and % N {iu) = 

\ V2ct n J 

Taylor theorem, for any |w| < 7r and sufficient large N, it holds 

6V2x(Jn (cot —Y- -_ _q^ 2 X(Jn f — w 


l/2 < 7 N 

-N 


ln/jv(w,ic) = — t , . — 

\ V 2 on 

Therefore, for N —> +oo, we have 




3\/2<j n 


■ By 


+ 0(cr N ) ) . 


lim f N (uj,z)=e 2 

N—too 


By lemma FT~T1 we have 


lim <f> N (iw) = lim 
N —)-+oo N—>oo 


sin 


—N 


V2(tn 


VZcn 
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Therefore By Theorem 14.11 it holds 


equivalently, 


l im I tt) P^{-2V0Nx) = 

N —>oo V iV / 777-! 


lim I 2 P»(-2VWx) = 

N—too \ iV / 777.! 


Corollary 4.2. 


lim M-f C» (J=) = 


N—too 


n—0 


Proof. By the generating function of Cff (x ), 

OO 

(1 - 2 xz + z 2 )- N = Y C^{x)z m , 

180 it holds 

OO / \ 

b (2ATlc ~ bl? 


Z = 


= i- 


l- 


XZ z 

7n + 2 N 

2 xz — z 2 
2 TV 


-N 


(4.4) 


□ 


(4.5) 


^(xz-z 2 / 2) 


Let = 


-1 / 2tc.z —z 2 \ 

1 l 2iV / 


, then limiv->-oo 9n = e. Therefore 


oo / \ 

S (2W, '" C ” U*) 


_ n xz n ~ z / 2 
* — 9 n 9 n 


and 


lim <&* = e xz 

N—too 

lim g,f' 2 = e- 1 / 2 . 


N—too 


By Theorem 4HD, we have 

^J2N)^C- (-|=) = ±H mM . 


□ 
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5. Biorthogonal systems relate to Generalized Laguerre polynomials 


Laguerre polynomials, L n (x), are solutions to the Laguerre differential equa¬ 
tion 


xy" + (1 — x)y' + ny = 0 , n > 0 . 


Laguerre polynomials is a class of orthogonal polynomials with weighting func¬ 
tion w(x) = e~ x . The Rodrigues representation for the Laguerre polynomials 


190 IS 


X J71 


and the generating function for Laguerre polynomials is 

OO 

(1 _ L m (x)z m , \z\ < 1. 


(5.1) 


(5.2) 


m —0 


The generalized Laguerre polynomials, Lm\ are also a class of orthogonal 
polynomials with weighting function w(x) = x a e~ x and generated by 


(1 - 2)-a-lgT=T = £ |z| < 1. 


(5.3) 


m =0 


The Rodrigues representation for the generalized Laguerre polynomials is 


rp—OLpX pjn 

-^(e-v+«). 


(5.4) 


195 When a = 0, we have L„\x) = L n (x). 
The explicit formula for L^\x) isQ 




™ (“ + !)i J ! 


(5.5) 


where (a) n = nr=i( a + * — 1), (a)o = 1, for n = 1, 2, 3_When a = 0, the 

explicit formula for generalized Laguerre polynomials, L^\x), become 


k =0 v 7 


(5.6) 


By simply computing, we have 
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Proposition 5.1. 


L 'n( x ) = L n-l(x) ~ L n -l(x). (5.7) 

200 The orthogonality relation for the Laguerre polynomials is contained in 

^ L^\x)L^\x)x a e- x dx = T{a + T ! + 1) S rnn , a > -1, (5.8) 

J o nl 

which can be considered as a biorthogonal relation between the derivatives of 
{-£k r(a+n+i) I and Laugerre polynomials { L:m = 0,1,...}, 


Corollary 5.1. 


T (“) 

1 


dx n T(a + n + 1) 


= Sr, 


lim {-l) m (2N)- m / 2 L<£\xV2N + N) = —H m {x). 
<f-ycx> ml 


Proof. 

°o / ~ \ ~ N ~ 1 "*7^ 

J2(-l) m (2N)- m/2 L£XxVm + N)z m = (l + -^=) e 1 + z /' /2N e 1+z ^' /2N . 

rrt.—O \ V / 


_ VNz 

~ za C: / \ iV+1 ^2 

/at(x,z) = e 1+2/V2IV , (f N {iz) = (1 + -7=1 gi-M/v^ then for iV — >• 00, we 
have 


linpv-i.00 /jv( a:, z) = e zx , 

2 2 

limAr-^oo <p N (iz) = . 

205 By Theoren i4.ll we have 

lim (-l) m (2N)- m / 2 L^\xVm + N) = —H m (x). 

N^roo ml 

□ 

In this section, we consider a family of biorthogonal polynomials, {P m (x, a, u>) : 
m = 0 , 1 ,...}, generated by a sequence of functions, <p(x, a , w), which converges 
to the generalized Laguerre polynomials, Lm\x). 
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Taking a compactly supported distribution <j> £ £'(K), Let (j) denote the 
Laplace transform of <f>. Then for any integer n > 0, 


(V n) , (1 - z) n e'-*^ = =z n '<i>{y^-) 


if m ^ o, 


-0 (1 - z) n e' : * \ = z „ 


(5.9) 


(5.10) 


\ *(£;) 

in a neighborhood of 0. Since <fi is compactly supported, <j> is analytic. So we 
can define a sequence of polynomials, P m , by the generating function 
(1 - z) n e ^ 

m =0 


= ^ P m(x)i 


(5.11) 


It follows from (15.101) and (15.111) that for any integer n > 0, 

OO 

z n = J2^ n \P m (x))z m , (5.12) 

m =0 

which gives the biorthogonal relation 

^ n \P m {x)) =6 m , n . (5.13) 

210 Differentiating (15.111) with respect to x and equating coefficients of z m in the 
resulting equation gives 

P' m ( x ) = P^x) - P m -i(x), m = 1,2 ,..., (5.14) 

which is similar to the property of Laguerre polynomials ^Proposition 15. ll) . 
^ J mi x ) = — L m —i(x). 

If 4>(x,a) = ^(a+n+i) > th e Laplace transform of<f)(x, a) is 4>(z,a) = ( z+1 )1+t»+i • 
which implies </>(yz-r,a) = (1 — z) a+n+1 . We have 

1 e~^(l -z) vX 


cj> {n \x,a), 


(5.15) 


0(r=I>a) / 

in a neighborhood of 0. We can define a sequence of polynomials, Lm\x ), by 
215 the generating function 

p— _ -y\ n _°° 

- ~r~ = (1 - z)~ a ~ 1 e~ T ^ 7 = J2 L^\x)z m . 

m=0 
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Q, —|_ ^ - U 

So the biorthogonal systems generated by function cf(x, a) = ^(a+n+i) are La- 
guerre polynomials. 


Theorem 5.1. Let (f>(x,a,uj) satisfy the following conditions: 

(1) There are constants 0 < r < 1 and c, for any e > 0, 36, such that 
220 |w — c| < S, it holds 


* rd- 


U - (1 - 2) 


Q!+n+1 


< e, 


2 < r. 


(5.16) 


Equivalently, lim u ^. c (j>(x,a,uj) = ■ 

(S) Lei {P m (:r, a, w) : m = 0,1,...} fee i/ie biorthogonal polynomials gener¬ 
ated by the functions, <j>(z, a, w) by 


e 1 ”*(l- 2 ) r 


= ^ Pm(x,a,uj)i 


1-z ’ a ' W ) m=0 

T/ien /or eac/i m = 0,1,..., {P m (x,a,w) : m = 0,1,...} converges locally 
uniformly to the generalized Laguerre polynomial, Lm\x ), as to goes to c. 

Proof Since (f>(0,a,uj) = 1, we can choose a neighborhood U of the origin so 
that > \ and |(1 — 2 )“ +1 | > | for all z £ U. Take a circle C 

completely contained in U, with centra at the origin 0 and radius r, so that 
(15.161) is satisfied. The coefficients of the Taylor series 
(1 — z) e x ~ z e 1 ~ z 


™=o 

are represented by the Cauchy’s integral formula: 

1 


E (Pm(x,a,L 0 )-Llz\x))z m (5.17) 


Pm(x,a,u) - L^\x) = 


Pm{x,a,uj) - L^\x) 


< 


< 


2i rz 

1 

2tt 

1 

2n 


-«(1 - Z ) a+n+1 - 
g 1 —2 -—-- QjZ 

C z m+1 <f>(j^,a,u))(l — z) a+1 

(l- Z ) a+n+ 1 -Hj^,a,u:) 


y.m-\- 1 


-dz 


1C r m+l 

< 4e xRe( -^£. 


K 1 — z) a+1 \ 
dz 


e xRe{ r ^) £ 


<P(r^,a,w) |(i-2)“ +1 | 
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It follows that for each m, P m (x,a,u>) converges locally uniformly to the gen¬ 
eralized Laguerre polynomial, Lm\x), as w goes to c. □ 


For the Meixner-Pollaczek polynomials, we have the generating function: 

OO 

F(x, z) = (1 - e^z) ~ X+ix (1 - e~ iu z) = £ P$> (*; u)z m . (5.18) 

n =0 


230 


Corollary 5.2. The Laguerre polynomials can be obtained from Meixner-Pllaczek 
polynomials by the substitution A = ^(a + l),x —> —^lo~ 1 x and letting u> —> 0. 

limP^(^; W )=L^(^). 

UJ—¥ 0 ZUJ 


Proof. 



Let ) = 


lim d>( 

w-K) 


z 


1 - 0 


[(1 - ze iu ) (1 - ze~ iu )] “= 1 (1 - z) n , it holds 
,a,w) = lim Ih-ze^il-ze-™)] 2 ^ (l-z) n 

= lim [l — 2z cos ui + z 2 l 2 (1 — z) n 

uj->0 L J 

= (1 - z) a+n+1 . 


Since 


lim 


/ 1 - ze iu> 
\ 1 — ze~ lul 


lim 

1 + z 

w-K) ' 

V 

lim 

'i + . 

w-K) ' 

V 

zx 

e 1 -*, 



e~ tu - e lw 
1 — ze~ iui 

e -iu _ e iu 

1 — ze~ lul 


1 — e %u> _ zx sin ui 
-2 iz sin u ' (l — Z e~ 


by theorem (E3D, we have 


lim P r 
)■() 


i+i -£ 
K 2uj 




L^\x). 


□ 
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(l - (1 - z)- p - x = E { -^M n {x-,(},c)z" 


235 The generating function for Meixner polynomials M n (x] /3,c) is 

(P)r. 

71—0 

Corollary 5.3. The Laguerre polynomials can be obtained from Meixner poly¬ 
nomials by the substitution ft = a + 1, x —> jff- and letting c —> 1 . 

^ 1 “' (a:) 


Proof. 


,, , cx 
Inn M u 

c-> l V 1 — c 


; a + 1, c = 


ii a) (0)' 


E 


71=0 




1- £\ — 


CX 

(1 — z ) —a— 1 

c ( l -*)-"- 1 


= 1 + 


1-z 

z(c- 1) 


z(l - z) 

(1 — z)c — zx 


Since 


(l-»)c 

.. , 1 z(c — 1)\ 1 ~ z _£_ a 

lim 1 + —-- = e 

c-s-l \ c(l - 2)7 


Therefore, by th.eorem lj5.lll . it holds 


lim(a+l ) n M n 
c-> i V 1 — c 


;a + 1 ,c ) = 


Since 


£W(x)= (a + 1)n E 


n! ft ( a + ’ 


2 « then £“(0) = (a + 1)„. Therefore 


.ex \ L^n\x) 

Inn M n - -; a + 1, c ) = 

c-n \ 1 — c 


4 a) (0)' 


(5.19) 


□ 
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